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In this work we give a consistent picture of the thermodynamic properties of bosons in the Mott 
insulating phase when loaded adiabatically into one-dimensional optical lattices. We find a crucial 
dependence of the temperature in the optical lattice on the doping level of the Mott insulator. 
In the undoped case, the temperature is of the order of the large onsite Hubbard interaction. In 
contrast, at a finite doping level the temperature jumps almost immediately to the order of the small 
hopping parameter. These two situations are investigated on the one hand by considering limiting 
cases like the atomic limit and the case of free fermions. On the other hand, they are examined 
using a quasi-particle conserving continuous unitary transformation extended by an approximate 
thermodynamics for hardcore particles. 
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I. INTRODUCTION 

The design of tunable quantum systems with strong 
correlations has attracted much interest due to the fast 
experimental developments in the field of ultracold atoms 
loaded in optical lattices. It has become possible in ex- 
periment to vary the dimension of the optical lattices 
and the particle statistics of the involved ultracold atoms 
QiHilE^IEIll' This opens the road to study fascinating 
quantum many-body effects like superfluidity in artifi- 
cially created and thus tunable structures. 

In order to investigate these quantum effects the tem- 
perature of the ultracold atoms in the optical lattice has 
to be small in units of typical model parameters. Exper- 
imentally it is an open problem to measure this tempera- 
ture. We have addressed this question in a recent work 
theoretically by analyzing spectroscopic experiments on 
bosons in one-dimensional optical lattices 3]. We found 
that the temperature in the Mott insulating phase is sur- 
prisingly large, namely of the order of the onsite Hubbard 
interaction U. In an analogous experiment the so-called 
Tonks-Girardeau regime (TO) was investigated ^ . From 
this experiment, a rather small temperature of the order 
of the hopping parameter t has been deduced. 

In the present work we provide a consistent thermo- 
dynamic picture of both experiments. The large differ- 
ence in temperature can be understood solely from en- 
tropy considerations. In the undoped Mott insulator the 
groundstate is unique and does not carry any entropy. It 
is separated from excited states by an energy gap which is 
determined by the interaction U . The entropy increases 
appreciably only at temperatures which are of the order 
of the gap. In a doped Mott insulator describing the TG 
regime the low-lying states produce a finite entropy and 
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the temperature of the bosons is much smaller. 

The paper is organized as follows. We will first intro- 
duce the model and our notation. In Sect. IIIII we will 
study several analytically solvable limits. First, we cal- 
culate the entropy of non-interacting bosons in a trap rel- 
evant for the cold atoms before loading into the optical 
lattice. Then we consider the atomic limit and the limit 
of non-interacting fermions which corresponds to the TG 
case. In Sect. IIVI a continuous unitary transformation 
(CUT) is used to study the Mott insulating phase of the 
Bose-Hubbard model at zero temperature. The results 
are extended to finite temperatures by an approximate 
formalism appropriate for hardcore particles. We first 
test the reliability of this approximation by comparison 
to exact results in the case of non-interacting fermions. 
Then the thermodynamic properties of the doped and 
undoped Mott insulating phase of the one-dimensional 
Bose-Hubbard model are studied. The conclusions will 
summarize the essence of the paper. In the last section, 
we propose an experimental procedure to reach very low 
temperatures also for the Mott insulating phase. 



II. MODEL 

Bosonic atoms in an optical lattice are described by 
the Bose-Hubbard model 1]. The Hamiltonian reads 

H = tHt + UHu + f-iHf, (1) 

i 

+ (C//2)^fi,(ni - 1) -M^fii 

i i 

where the first term is the kinetic part tHt and the second 
term the repulsive interaction UHjj with U > 0. The last 
term accounts for a chemical potential /i. The bosonic 
annihilation (creation) operators are denoted by b^^\ the 
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number of bosons by fii = The average number of 

bosons is n = {fii) = N/L where N is the total number 
of bosons and L the total number of wells in the optical 
lattice. In the spirit of the tight-binding description we 
will call the wells henceforth 'sites'. On the same site, 
the bosons repel each other with interaction strength U. 

In the limit t/U 0, a site with n bosons contributes 
the energy E{n) = {U/2)n{n — 1). A deviation from 
uniform filling costs the energy E{n + 1) + E{n — 1) — 
2E{n) = U . Note that this excess energy results from 
the constant curvature of the parabola oc n^. Thus it is 
independent of n. For large U and integer filling n it is 
energetically favorable to distribute the bosons uniformly 
and to put just n bosons on each site. The system is in 
an insulating state with an energy gap of the order of U . 

In the opposite limit of large t/U the bosons form a 
superfluid condensate for low enough temperatures. 

We will not address the inhomogeneity that is present 
in the experimentally realized traps, i.e., we treat the sys- 
tem as translationally invariant though there is the slowly 
varying external trap potential. This approach makes 
quantitative computations much easier and it is backed 
by the recent finding that the inhomogeneous system acts 
essentially like two (or several) independent systems of 
different average number of bosons [g] ■ 

Among the extensive literature on the bosonic Hub- 
bard model there are mean-field treatments, e.g. Refs.Q, 
11, ^pl. other analytical investigations, e.g. Refs. 
1 3. [ iSl ll4L fill and numerical treatments, e.g. Refs. 
7l lid llT l li a Il9i The phase diagram in the t-^-plane 
consists of a series of lobes. The lobes of the Mott in- 
sulating phases correspond to an integer filling per site. 
The present paper is only concerned with the first lobe 
of the phase diagram where n = \. 



III. LIMITING CASES 



Three limiting cases are investigated in order to un- 
derstand the basic entropy properties of the cold atoms 
before and after adiabatic loading into the optical lattice. 
First, we calculate the entropy of non-interacting bosons 
in a trap. This entropy is present before turning on the 
optical lattice. It can be viewed as a function of the su- 
perfluid fraction. The same entropy has to be realized 
in the optical lattice once it has been turned on in an 
adiabatic fashion. Hence the adiabaticity of the loading 
determines the temperature of the bosons in the optical 
lattice [23. We calculate the entropy as a function of 
the microscopic model parameters and of the doping in 
the atomic limit, t/U = and in the limit of infinite U , 
the so-called TG regime. In one dimension I21I |2^ [23l|. 
a Jordan- Wigner mapping to non-interacting fermions is 
appropriate. 



A. Free Bosons 

The formula for the entropy of free bosons in a wide 
harmonic trap is derived. The Hamiltonian of the har- 
monic oscillator in d dimensions reads [h — \) 



(2) 



Many bosons loaded in the harmonic trap are considered. 
The energy ujq is small. Therefore, sums over multiples 
of loq can be rewritten as integrals. The groundstate is 
macroscopically occupied with a fraction 



/o = No/N 



(3) 



where N is the total number of bosons and A^o their num- 
ber in the condensate. With the density of states D{lu) 
and the boson occupation riB(w) = l/(e'^" — 1) one has 
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nB{mujo)D(mLL}Q) (4a) 



nBiuj)D{Lu)duj . (4b) 



The density of states D{llj) counts the number of states 
in a given energy interval. Let us consider intervals of 
length loq centered around the energies muQ, where m is 
a non-negative integer. We consider the number of ways 
in which the energy u = rriLOQ can be realized. One has 
to distribute m indistinguishable energy quanta over d 
oscillators. The number of possibilities is 



i + d-l 
d-1 



Using m = uj/ivo, the leading power in lo reads 



D{u;) = 



1 



LOoid - 1)! \ujo 



(5) 



(6) 



Only the leading order in m — lu/loo is retained because 
a wide trap potential is considered which implies a small 
oscillator frequency ujo- 

Specializing to our case of interest, a three-dimensional 
trap (d = 3), Eg. I4bl becomes 



N^Nfo + 



1 



dx = Nfo+T^C{3), (7) 



where x = (3u! and C(z) is the Riemann ('-function. Fur- 
thermore, we need the total entropy. The entropy of a 
single bosonic level e reads (fee — 1) 



S{e) = -\n{l - e-^') + PeuBie) 
The total entropy is given by 

5- = r S{e)Die)de = ^ f- 
Jo 45 V^o 



(8) 



(9) 
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The condensate at zero energy does not contribute to 
the entropy because it represents a single, pure quantum 
state. Combining Eqs. and 10 leads to 



S/N = (1 - /o) 



2^4 



45C(3) 



3.60(1 - /o) 



(10) 



This relation between entropy and superfluid fraction will 
serve in the following as the initial entropy of the cold 
atoms before turning on the optical lattice. 



B. Atomic Limit 

In this part, we calculate the entropy as a function 
of temperature and doping in the atomic limit, namely 
t/U — 0. Then, all excitation processes are completely 
local. For simplicity, we restrict the local Hilbert space 
to the state with occupation one and the two energeti- 
cally adjacent states with occupation zero and two. The 
energy to create a site with two bosons is U. All states 
with an even larger number of bosons have a higher en- 
ergy and are therefore not important for the entropy in 
the temperature range considered here. 

The Gibbs energy of a single site is given by 



= -rin 1 + e/'A' + e-'3(t^-2,0 



(11) 



where the terms in the argument of the logarithm refer to 
the empty, the singly occupied and the doubly occupied 
site. Introducing the doping S := n—1 and using dfj_g\rp — 
—n = —[1 + 5) we arrive at the following expression for 
the chemical potential at given doping 5 



2(1-5) 



(12) 

The entropy per particle is calculated as the derivative 
of the Gibbs energy with respect to temperature S/N — 
S/{nL) = — drfl^- In Fig. ^ the entropy is shown for 
n e {1; 0.99; 0.95; 0.9; 0.8}. Additionally, the entropy for 
free bosons with a superfluid fraction /o = {0.95; 0.9; 0.8} 
according to Eq.EHis given. All curves display an almost 
constant plateau at small temperatures T/U < 0.1, an 
uprise at about T/U 0.3 and a saturation for large 
temperatures. 

Supposing that the bosons are loaded adiabatically, 
three different regimes can be identified for fixed /q. In 
the following, the experimentally relevant value /o — 0.9 
is discussed. At zero doping, the temperature is ap- 
proximately U/6, i.e., of the same order as deduced by 
analysing recent spectroscopic experiments 3, 2]- This 
seemingly large value for the temperature is needed be- 
cause the groundstate of the Mott insulating phase at 
zero doping is separated by a large gap from the first 
excited states. The entropy has to be produced by the 
excited states which have an energy of the order U. 
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FIG. 1: Entropy per particle as a function of temperature in 
the atomic limit t/U = for various values of doping. The 
entropy per particle of the three-dimensional harmonic oscil- 
lator for condensate fractions /o = 0.95; 0.9; 0.8 is included as 
black horizontal lines. 



A small amount of doping leaves this temperature al- 
most unchanged. But the system develops a finite en- 
tropy at zero (and small) temperatures. The reason for 
this is the degeneracy of the groundstate coming from the 
different possibilities to distribute N bosons over L lat- 
tice sites. The entropy at intermediate values of the tem- 
perature is only changed slightly. The temperature real- 
ized in the optical lattice is changed abruptly when the 
zero temperature entropy reaches the entropy value of the 
free bosons. Then the temperature immediately jumps to 
zero. Physically this corresponds to the situation where 
the whole entropy can be produced by the groundstate 
degeneracy of the doped Mott insulator. Thereby the 
temperature is reduced to the scale t which is zero in 
the atomic limit. This is the mechanism which is most 
probably observed in the experiment in Ref . 0. It is more 
closely examined in the next part. 

In the third temperature regime, the entropy of the 
free bosons is smaller than the entropy of the adiabati- 
cally loaded bosons at any temperature. This means that 
an adiabatic loading of the bosons into the optical lat- 
tice is not possible. Instead, a non-ergodic mixture will 
be realized. We emphasize, however, that this feature 
results from the macroscopic ground-state degeneracy in 
the atomic limit which represents an extreme and ideal- 
ized situation. 



C. Tonks-Girardeau Regime 

The experiment performed in Ref. |^ probes the so- 
called Tonks-Girardeau (TG) regime HI H H. 
This is the limit of large U. Only empty and singly oc- 
cupied sites are present and the problem can be mapped 
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FIG. 2: Entropy of non-interacting tight binding fermions 
corresponding to the TG regime for various fillings. The en- 
tropy per particle of the three-dimensional harmonic oscillator 
for condensate fraction /o = 0.95; 0.9; 0.8 is included as black 
horizontal lines. 



in one dimension onto fermions without interaction if we 
consider nearest-neighbor hopping only. If we relax the 
condition of nearest-neighbor hopping the mapping to 
fermions can still be done but the resulting system is 
not without interaction. Hence we focus here on nearest- 
neighbor hopping so that we have a benchmark system 
which is simple to analyze. 

In order to know which temperature is realized after 
adiabatic loading in the TG regime one needs to know 
the entropy of the non-interacting fermion gas. In one 
dimension the dispersion of tight-binding fermions with 
hopping ~t reads 



e(fc) = -2tcos(fc). 
The Gibbs energy per site is (k^ = 1) 

G 
L 



5 = T 



(13) 



(14) 



(15) 



where L is the number of sites. The entropy per particle 
is calculated as 



N 



IS 
n L 
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IT 



1 dg_ 

n df 

ln(l4 



(16) 



-/3e(fc) 



) + 



6(fc) 



dk 



T{l + e-rMk))^ 

In the limit of high temperatures the entropy per site is 

S/L^-n\og{n)-{l~n)log{l-n) (17) 

which is just the entropy of one state with occupation 
probability n. The entropy per particle is accordingly 



S/N = - log(n) -I- (1 



^)log(l-n). (18) 



Figure 12 shows the entropy per particle versus T/t in the 
TG regime for fillings n € {0.5,0.6,0.7,0.8,0.9}. After a 
steep increase the entropy reaches its saturation value of 
— log(n) + (1 — n~^) log(l — n). Filling n = 0.5 shows the 
largest saturation value. If the filling is increased further 
the saturation value decreases. Note that the S/L is 
particle- hole symmetric, i.e., it does not matter whether 
we consider the filling of fermions n or the filling of holes 
1 - n. 

Additionally, the entropy per particle of the three- 
dimensional harmonic trap is plotted for condensate frac- 
tions /o = 0.95; 0.9; 0.8 as horizontal lines. This allows 
to determine the temperature that results for these con- 
densate fractions after loading into the optical lattice in 
the TG regime. For half filling n ~ 0.5, the temper- 
atures are T/t = 0.17; 0.32; 0.61 for /o = 0.95; 0.9; 0.8. 
Also temperatures of multiples of the hopping parameter 
are possible, e.g., for /g = 0.838 and n = 0.8 the tem- 
perature is T/t = 2. Again, the temperatures are of the 
size of the microscopic system parameter t. This agrees 
well with the temperatures estimated in Ref. 4. There, 
the temperature has been determined from the measured 
momentum profiles. 

For higher filling the entropy saturates at a lower value. 
For filling n = 0.8 it does not rise up to the value of the 
harmonic trap for /g = 0.8. Thus in a process of adia- 
batic loading one cannot achieve a filling of n = 0.8 when 
starting with a condensate fraction of /o = 0.8. Within 
the TG regime, even the highest temperatures are not 
sufficient to reach the necessary entropy. In practice this 
means that a system starting from too high an initial 
entropy will not end up in the TG regime. Doubly occu- 
pied sites will occur corresponding to temperatures of the 
order of the corresponding gap. This means that the de- 
scription as non-interacting gas of fermions will no longer 
hold. 



IV. CONTINUOUS UNITARY 
TRANSFORMATION 



After having discussed the atomic limit and the case 
of free fermions, we study the Bose-Hubbard Hamilto- 
nian at finite t/U by means of a particle conserving CUT 
at zero temperature 0. We reformulate the model in 
terms of new quasi-particles above the reference state 
with one boson per site. These particles obey hardcore 
statistics. Aspects of the method will be presented in 
the next subsection. The focus is again laid on ther- 
modynamic, mainly entropy, considerations. We adapt 
a formalism introduced for spin ladders psj to approx- 
imate the hardcore statistics to obtain results at finite 
temperature In subsection IIV B| this formalism is in- 
troduced for the case of the Bose-Hubbard Hamiltonian. 
In particular, we check the reliability of the approximate 
hardcore statistics by comparison to the results of non- 
interacting fermions in the TG regime. At the end of this 
part, we present results for the entropy and the temper- 
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ature of the undoped and doped Mott insulating phase 
for finite t/U . 



A. Method 

In this part some technical aspects concerning the CUT 
are presented. For the general method we refer the in- 
terested reader to Refs. H IHilii To set up the 
CUT calculation the reference state has to be fixed. For 
t/U and filling n = 1 the groundstate of Eq. ^is the 
product state of precisely one boson per site 



|ref) 



N, 



(19) 



where \n)i denotes the local state at site i with n bosons. 
We take |ref) as the reference state. For finite t/U the 
kinetic part tHt will cause fluctuations around the ref- 
erence state. Deviations from |ref) are considered to be 
elementary excitations. Each site can be occupied by an 
infinite number of bosons. The local Hilbert space has 
infinite dimensions. Therefore, there are also infinitely 
many linear independent local operators on this space. 
To set up a real space CUT, it is convenient to truncate 
the otherwise proliferating number of terms 0, ll^l ■ In 
the present work, the local bosonic Hilbert space is trun- 
cated to four states. Numerical studies using density- 
matrix renormalization group have shown that this does 
not change the relevant physics The CUT 

gives reliable results in a wide parameter regime up to 
t/U ~ 0.2, see Ref. J,. Since the calculation is restricted 
to four local states including the reference state there 
are three local creation operators. Applied to the refer- 
ence state |ref)i = they create local excitations. A 
hole at site i is created by = \0)i, a particle by 

pl\l)i = \2)i, and = |3)i induces another kind of 

particle at site i. The particle created by dj corresponds 
to a local occupation of a site by three bosons which is 
treated as an independent quasi-particle. The operators 
h, p and d obey the commutation relations for hardcore 
bosons. The local bosonic operator is thus represented 
by the 4x4 matrix 



V2pl 



/ \ 

10 

\/2 

V V3 / 



(20) 



which acts on the vector of the coefficients of 
{|0), |1), |2), |3)}. The CUT is implemented using this 
presentation of the Hamiltonian with finite local matri- 
ces [23 • The Hamiltonian Eq.^ conserves the number of 
bosons. But if it is rewritten in terms of /i, p, d it does not 
conserve the number of these newly introduced particles. 
For example, the application of Ht to |ref) generates a 
particle-hole pair. To eliminate the terms that change 
the eigenvalue of Hu an infinitesimal generator of the 
CUT is chosen as 



in an eigenbasis of Hu; q,; is the corresponding eigen- 
value of Hu, see Ref. |2^ The choice of this generator 
leads to an effective Hamiltonian commuting with Hu 
which implies that the processes changing the number of 
quasiparticles are eliminated. For the bosonic Hubbard 
model a complete elimination cannot be achieved due to 
the occurrence of life time effects Q • But this effect is of 
minor importance in the parameter regime under study. 

The bottom line is that the CUT maps the complicated 
elementary excitations of the bosonic Mott insulator onto 
hardcore particles. The corresponding annihilation (cre- 
ation) operators are p^^' , h<-^\ and d(t). 



B. Finite Temperatures 

In order to discuss thermodynamic properties of the 
Bose-Hubbard model we have to extend the calculation to 
finite temperatures. The CUT yields an effective Hamil- 
tonian. The particles involved obey hardcore commuta- 
tion relations. The statistics of hardcore particles, how- 
ever, is complicated and thus a direct evaluation of the 
partition function is not easily tractable. In order to ob- 
tain results at finite temperature we will adapt the cal- 
culation in Ref. 25 for hardcore triplet excitations, called 
triplons, to our present purposes. 

The explicit derivation of the approximate Gibbs en- 
ergy for hardcore bosons is presented in Appendix ^ 
The resulting Gibbs energy per particle reads 



g = -{T/L) \nZ = -T ln(l + zP + z'' 



(22) 



where we use partition functions for the individual par- 
ticles defined in Eq. (|A5|) . 



Vi,3 (0 = sgn {qi - q.j) Hi^j (l) 



(21) 



C. Test of the Approximate Thermodynamics 

We check the validity of the approximate thermody- 
namics. The results for the TG regime from the previous 
section. Sect. IIII CI are used as a benchmark. In par- 
ticular, the quality of the approximate statistics for low 
temperatures and fillings n 7^ 1 is addressed. 

For this purpose, we introduce a toy model which cor- 
responds to the TG regime. The CUT calculation deals 
with the situation found for the bosonic Hubbard model 
for filling n = 1 and for values t/U in the Mott insu- 
lating phase. Above the groundstate of singly occupied 
sites there exist particle- and hole-like excitations. The 
toy model includes only the hole-like particles. Its Gibbs 
energy in Eg. 1221 simplifies accordingly (cf. Appendix A) 



5toy = -Tln(l + z'') . (23) 
The dispersion contains only a single cosine-term 

uj^{k) ^ -2tcos{k) (24) 
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FIG. 3: Comparison of entropies per particle at small tem- 
peratures calculated for the TG regime, (a) (color online) 
Entropy of free fermions is shown as dotted black curves. Re- 
sults from the toy model calculation in gray (cyan). For fill- 
ings close to n = 1 the entropies show good agreement. But 
deviations increase with increasing doping and for decreasing 
temperature, (b) Relative deviation As = {Stoj — Stg)/Stg 
of the entropy of the toy-model from the TG result. 
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FIG. 4: (color online) The toy model calculation differs from 
the exact fermionic calculation at low temperatures. The tem- 
peratures T'^^^ and characterizing the limit of the ap- 
proximate thermodynamics are plotted versus = [1 — n)^ . 
For comparison, T"^''"' resulting from the thermal de Broglie 
wave length A is included. For further details see main text. 



behavior is still found for very low T. 

To investigate systematically, below which tempera- 
ture the approximate thermodynamics becomes signifi- 
cantly wrong, we define the temperature T'^°'^ , at which 
the relative deviation AS/S exceeds 10%. It is displayed 
in Fig. ^ as a function of 6^ = {1 — n)^. In addition, the 
temperature r"°s jg shown which is defined as the tem- 
perature where the entropy of the toy model becomes 
negative. The temperatures T"^^^ and T'^°s characterize 
the limit of applicability of the approximate thermody- 
namics. 

Both temperatures are compared in Fig. ^ to the tem- 
perature T^^'^^ which results from equating the thermal 
de Broglie wave length 



corresponding to hardcore bosons with nearest-neighbor 
hopping only. With these definitions the toy model de- 
scribes the TG regime. 

The results for this toy model calculation are depicted 
in Fig. Ell in comparison to the entropy from the exact 
formula for non-interacting fermions. Discrepancies are 
found to be increasing with doping. They appear espe- 
cially at low temperatures. 

Since the deviations are small the relative deviation 
AS/S = (5toy - 5'tg)/5'tg is plotted in Fig. EJs. For 
n ~ 0.8 the difference of the entropies has a maximum at 
T/t 0.08, then it drops and becomes even negative at 
very low temperatures because Stay becomes negative for 
T — > 0. This is an unphysical artifact of the approximate 
thermodynamics for hardcore bosons. For n ^ \ the 
maximum of the relative deviation is lower and shifted 
to lower temperatures. But artificial negative entropy 



with the average distance 1/5 between two hole excita- 
tions h. The inverse mass l/m is given by 2t which is 
the value of the second derivative at the band minimum. 
The excellent agreement in the proportionality T oc 6^ 
and the very good agreement of T'^^'^^ and T'^°^ clearly 
show that the approximate thermodynamics is a classi- 
cal treatment, which is valid as long as the gas of excita- 
tions is dilute enough. This argument quantifies the lim- 
itations of the approximate thermodynamics. For very 
low temperatures the results become unphysical. On the 
other hand, it shows that the results obtained for inter- 
mediate and high temperatures agree very well with the 
exact behavior in the TG regime. 
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D. Results 

In Sect, mil we have already learned that the tem- 
perature of the bosons in the optical lattice depends 
strongly on the regime of the system (Mott insulator or 
TG regime). The crucial point is whether the energeti- 
cally low-lying states realize the initially present entropy 
or not. In the following, we investigate the entropy at 
finite t/U for the effective model obtained by the CUT 
calculation and evaluated at finite temperatures by the 
approximate statistics for hardcore bosons. 

In the limit of small t/U , the CUT can be performed 
straightforwardly. The results for the Mott phase are ob- 
tained for commensurate filling n = 1. For fillings n ^ 1 
the reference state used for normal ordering before trun- 
cation should ideally be the disordered state at the lower 
filling 71 7, 26J. This requires to redo the whole CUT 
for each n. But the effect of a different reference state 
is not very large if t/U is small. So our approach is to 
perform the CUT with the reference state (I19|l at n = 1 
as described in detail in Ref. 0- Note that this approach 
is exact as long as the CUT is performed without any 
approximation. So it does not represent an additional 
approximation in itself. The deviation from rt = 1 is 
included only on the level of the approximate thermody- 
namics of hardcore bosons described in Sect. IIVBI By 
tuning the chemical potential in Ea. lASI the average num- 
ber of holes can be changed. This way of accounting for 
doping works particularly well for n w 1. 

The entropy per site for the Bose-Hubbard system is 
calculated from the Gibbs energy (|22|) a,s s — S/N — 

— ^ . First, the entropy in the TG regime is discussed. 

In Fig. loj the results for the entropy per particle from 
the CUT are compared to the results for the TG regime. 
The figure shows results for t/U = 0.01 and t/U = 0.02 
corresponding to large values of the parameter 7 = U/t 
where the TG gas is realized H The thick dotted 
curve shows the entropy per particle for t/U = 0.01 for 
filling n = 1. The increase of entropy with tempera- 
ture is much slower for n = 1 than for fillings n = 0.9 
(black curves) and n = 0.8 (gray/cyan curves). For filling 
n = 0.9, a very good agreement is found for low tempera- 
tures. The results based on CUT and on the approximate 
thermodynamics reproduce nicely the entropy in the TG 
regime. Good agreement is found up to temperatures of 
T/t « 10. For these temperatures, the TG result has 
already saturated at its high temperature limit. On the 
other hand, the CUT results begin to feel the presence of 
particle states (double occupancies in terms of the origi- 
nal bosons) which are not included in the TG description. 
Thus the CUT results deviate from the TG ones. The 
entropy is plotted versus T/t in Fig. |S1 A temperature 
of mt corresponds to T/U = 0.1 for t/U = 0.01 but 
to T/U = 0.2 for t/U = 0.02. This is why the resuhs 
for t/U = 0.02 show deviations already for T/t « 5. A 
similar behavior is seen for n = 0.8. 

The small deviations discernible at low temperatures 
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FIG. 5: (color online) Entropy per particle in the TG regime. 
For small temperatures, the results from the CUT and the 
free Fermion calculation agree. At higher temperatures there 
are deviations caused by the presence of double occupancies 
in the bosonic Hubbard model. 



(T « t) must be attributed to the difference between the 
effective model obtained from the CUT and the simple 
fermionic tight-binding model with nearest-neigbor hop- 
ping only in the TG regime. 

Finally, Fig. displays the entropy from the CUT cal- 
culation as function of temperature in units of U . The 
entropy per site for the undoped case is shown in Fig. 
|S1 for three values of t/U (black curves). The entropy 
does not change significantly for the values of t/U that 
are shown because all of them are still deep in the Mott 
insulating phase. 

The entropy per particle for condensate fractions /o = 
{0.95; 0.9; 0.8} is included as horizontal hues in Fig. 
Assuming such condensate fractions and a hopping pa- 
rameter t/U = 0.02 the corresponding temperatures are 
0.12; 0.16; 0.27 in units of the interaction parameter U. 
These values are within the range of temperatures that 
were found by the analysis of spectral weights 0. 

The inffuence of doping is shown in Fig. by cyan 
(gray) curves for n = 0.9. Clearly, several differences 
can be identified. At finite doping and small t/U, the 
constant entropy plateau seen for t/U = in Fig. ^ can 
still be recognized but there are changes for small tem- 
peratures of the order t. There is a finite dispersion at 
any finite t. Therefore, the entropy at zero temperature 
is zero because the macroscopic groundstate degeneracy 
is lifted. For increasing t the plateau is washed out. At 
very low temperatures, there is a sharp uprise from zero 
entropy to a reduced entropy value comparable to the 
value realized at the plateau in the atomic limit. At 
higher temperatures, the entropy is dominated by the 
particle states (double occupancies) in analogy to what 
was found in the undoped case. 
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FIG. 6: (color online) Entropy per particle as function of tem- 
perature for tjU — 0.02, 0.05, 0.1 at zero doping (black curves) 
and finite doping n = 0.9 (cyan/gray curves). In addition, the 
value for the entropy per particle of the three dimensional har- 
monic oscillator with condensate fraction fo — 0.8, 0.9, 0.95 is 
given. 



V. CONCLUSIONS 

The main purpose of this work was to clarify the differ- 
ent orders of magnitude of temperature in recent exper- 
iments on the Mott insulating phase of ultracold bosons 
loaded adiabatically in one-dimensional optical lattices. 
On the one hand, temperatures of the order U are found 
for the undoped Mott insulating phase [3, 13 • the 
other hand, temperatures of the order t are found in the 
TG regime. 

By thermodynamic considerations, we found a physi- 
cally relative simple and consistent picture. In the un- 
doped Mott insulator, the temperature has to be of the 
order U. The groundstate is a single state and does not 
provide any entropy. The initially present entropy has 
to be realized by excited states. This requires a large 
temperature of the order of the gap. The situation is 
different for the doped Mott insulator. Already at the 
energy scale of the hopping an appreciable amount of 
entropy is possible due to the finite amount of holes or 
particles. 

In comparison to experiment, an item still missing in 
our treatment is the inhomogeneity of the trap. So one 
might wonder whether our qualitative conclusions apply 
to the realistic systems. This question can be answered 
positively because it has been shown previously that the 
groundstate displays constant density plateaus at small 
values of t/U [SfljIsJl- For zero or low temperature, the 
realistic situation can be seen as an ensemble of homoge- 
neous systems of various fillings [||. Thus it is expected 
that our qualitative conclusions concerning the orders of 
magnitude of temperature remain valid even if the inho- 
mogeneity of a wide trap is taken into account. 

The results were obtained on the one hand by consid- 



ering limiting cases where a purely analytic treatment 
is possible and on the other hand by using a particle 
conserving CUT at zero temperature extended by an ap- 
proximate thermodynamics for hardcore bosons. For the 
CUT the Bose-Hubbard Hamiltonian is reformulated in 
terms of a reference state with exactly one boson per 
site and excitations above this state. This approach 
yields reliable results in a wide parameter regime of the 
Mott insulating phase. Results at finite doping were ob- 
tained by tuning the chemical potential in the approxi- 
mate hardcore thermodynamics. The reference state was 
not adapted to finite doping, which may represent a way 
to improve the CUT calculation. But we do not expect 
significant changes. 

We tested the reliability of the approximate thermody- 
namics for hardcore bosons by comparison to the exact 
results in the TG regime which are based on the map- 
ping of the system to non-interacting fermions. We find 
that the approximation works well up to relatively low 
temperatures T > t/2. We found that the approximate 
hardcore statistics corresponds to a classical description 
of the gas of hardcore particles. This results from the 
observation that the temperature below which the ap- 
proximate statistics fails can be estimated to be the tem- 
perature at which the thermal de Broglie wave length At 
equals the interparticle distance. 

Our considerations were restricted to one dimension 
because only in one dimension the hardcore bosons can 
be exactly mapped to fermions providing an easy-to-test 
benchmark. But apart from this generic one-dimensional 
feature used for the benchmark none of our arguments 
is restricted to one dimension. Thus we expect that the 
qualitative findings presented in this paper apply also to 
higher dimensions. 

The obtained results provide a consistent picture of the 
role of temperature for bosons in one-dimensional optical 
lattices. We found that different physical situations can 
produce rather different temperatures. It is important 
in future experimental and theoretical works to take the 
influence of temperature into account in order to explore 
new quantum phenomena in the field of ultracold atoms. 



VI. OUTLOOK 

Our basic finding is that a restricted phase space at 
low energies leads inevitably to a relatively large tem- 
perature when the system is reached adiabatically from 
a Bose-Einstein condensate with a condensate fraction 
which deviates from unity. One obvious means to reach 
lower temperatures is to realize condensates with less en- 
tropy, i.e., with a higher condensate fraction /o at lower 
temperature. 

Let us, however, suppose that this route cannot be used 
to reach the desired very low temperatures which are, for 
instance, necessary to perform reliable quantum compu- 
tations. In such a case, we like to propose an alternative 
route. The system with very restricted phase space at 
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low energies should be brought into weak thermal con- 
tact with another system in the optical trap which has a 
large phase space at low energies. Elementary thermody- 
namics tells us that the temperatures of both systems will 
finally become equal. This means that the system with 
large phase space cools the wanted system with small 
phase space down to low temperatures. 

In the realm of bosons loaded in optical lattices this 
proposal means that one system should be realized which 
corresponds to a Mott insulator and which remains at 
relatively high temperatures in the adiabatic process be- 
cause it does not have phase space below its gap. This 
Mott insulator can be cooled by a system which is in the 
TG regime, i.e., in the doped Mott insulating regime, 
with small hopping amplitude t. The TG system can ac- 
commodate for an appreciable amount of entropy even 
at low temperatures. Thus it is well-suited to cool the 
other system. 

We cannot devise an experimental system which real- 
izes the above idea. The use of two different laser systems 
in spatial vicinity for the same kind of atoms is certainly 
cumbersome. So it might be more promising to use two 
kinds of atoms in the same optical trap. The control pa- 
rameters for the first kind of atoms should be such that 
it becomes Mott insulating upon loading into the opti- 
cal lattice. The second kind should enter the TG regime 
upon loading. It is most challenging to see to how low 
temperatures a bosonic Mott insulating system can be 
pushed in such a way. 
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APPENDIX A: APPROXIMATE 
THERMODYNAMICS FOR HARDCORE 
BOSONS 

This appendix provides the derivation of the approx- 
imate Gibbs energy of hardcore particles used in Sect. 
IIVI The calculation adapts the reasoning in Ref. for 
triplets in a spin ladder to the present hardcore excita- 
tions. The relevant excitations are independent particle 
and independent hole excitations. Multiple correlated p 
and h excitations and excitations of d-type are less rel- 
evant because of their higher energy. Therefore, we will 
neglect d-excitations and correlated excitations of more 
than one particle or hole in the following. 

To use the non-interacting boson statistics for hardcore 
particles is a bad approximation because non-interacting 
bosons have no restriction on the occupation number. An 
incorrect number of states would contribute in the par- 
tition sum. The idea is to reweight all contributions to 



the partition sum such that they contribute the correct 
amount of entropy. Therefore we have to calculate the 
correct number of states with particle and hole excita- 
tions. 

The number of possibilities to distribute M truly 
bosonic excitations on 2L sites is 



9b{L,M) 



2L 



-M 
M 



1 



(Al) 



But the correct number of possibilities of M hardcore 
excitations of two kinds reads 



5(L,M) 



M 



(A2) 



To obtain the correct entropy for each number of excita- 
tions the Af -boson part in the partition sum is reweighted 
by the factor g{L, M)/gB{L, M). By this means, the par- 
tition function displays the correct high and low tempera- 
ture behavior. The boson partition function for a system 
with L sites and M bosons is 



Z'(L,M)= 



M=Mp + Mh 

(A3) 

The sum runs over all multi-indices {kj,ki} of Af = 
Mp + Mh momenta for particles at the momenta {kj} 
and holes at the momenta {ki}. Here, u)^{kj) is the en- 
ergy of a particle with momentum fcj, ui^{ki) the energy 
of a hole with momentum ki. The chemical potential /i 
is accounted for by /i^ = /i and fjL^ = — /i. The nega- 
tive sign in the definition of the chemical potential fi^ 
for holes results from the fact that adding a hole means 
removing an original boson. The boson partition func- 
tion Eg. IA3I contains (2L)^ terms. Thus the reweighted 
partition function reads 



A/=0 
L 



(A4a) 



M=0 



{kj.ki) 

M=Mp+Mh 



L 



/3((^-(fc)-p<') 



where we used the definition 



1 



27r 



27r 



(A4c) 
(A4d) 

(A5) 



for the partition function of the single modes with a G 
{p, h\. The chemical potential /i is determined by speci- 
fying a certain filling. For instance, for the filling n = 1 
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the chemical potential ^ is determined such that as many 
particles as holes are excited. 

We emphasize that Eq. IA4I is only an approximation 
to the correct hardcore partition function because it con- 
tains contributions from states that are only allowed for 
non-interacting bosons, but are forbidden for hardcore 
bosons. The correct contribution of each subspace of 
fixed total number M of excitations is ensured only in an 
overall manner by reweighting the total subspace. 

Using the final expression (|A4d|) one can compute the 
Gibbs energy per site g to be 
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- {T/L) In Z ^ -T\tl{1 + zP + 



(A6) 



Eauations lA5l and lA6l are the formulae used in Sect. II V 51 
The toy model discussed in Sect. IIV CI includes only 
holes as possible excitations. Thus the Gibbs energy de- 
pends only on one partition function of a single mode 



.9toy = — 'rin(l -I- z^). The dispersion of the toy model 
contains one single cosine-term uj^{k) = —2t cos(fc) which 
makes the analytical calculation of z'' possible 



(A7) 



where is the modified Bessel function of the first 

kind. The entropy per particle reads 



S_ 

N 



1 % 
n dT 



2th m) + M^^^o (2/3t) 

+ ilog(l + e'5^'/o(2/3i) 
n V 



(A8) 



(A9) 
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